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ABSTRACT 

A lattice definition of QCD based on chiral defect fermions is discussed 
in detail. This formulation involves (infinitely) many heavy regulator 
fields, realized through the introduction of an unphysical fifth dimension. 
It is proved that non-singlet axial symmetries become exact in the limit 
of an infinite fifth dimension, and before the continuum limit is taken. 
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1. Introduction 



Consistent local regularization methods which preserve gauge invariance must 
break axial symmetries explicitly. This is a consequence of the well-known axial 
anomaly |jl|]. However, flavour Non-Singlet Axial Symmetries (NSAS for short) are 
recovered in renormalized correlation functions to all orders in perturbation theory 0, 
^. (For recent progress and references to earlier literature see ref. Q). 

Going beyond perturbation theory, the rigorous definition of QCD relies on the 
lattice regulator. The most popular method to avoid the fermion doubling prob- 
lem employs Wilson's prescription for the fermion action p. The advantage of this 
method is in its simplicity. But the ensuing breaking of axial symmetries is hard, in 
the sense that perturbative corrections to quark masses are 0(l/a) where a is the 
lattice spacing. Hence, one has to fine tune the bare quark masses in order to recover 
the correct renormalized masses in the continuum limit. 

Using Wilson fermions, it was shown that weak coupling perturbation theory 
(WCPT) on the lattice reproduces the axial anomaly |0] , and that NSAS are recovered 
to all orders in WCPT in the continuum limit 0. These results have in fact some 
validity beyond the scope of WCPT, and one can discuss the renormalization of 
gauge invariant composite operators. But, because of the severe fine tuning problem 
inherent to Wilson fermions, one cannot give a completely general non-perturbative 
proof of the restoration of NSAS. 

Assuming the existence of the chiral limit in the full quantum theory, what one can 
do in numerical simulations with Wilson fermions is to determine the correct finely 
tuned values of the bare masses by measuring some correlation functions. Fixing 
the bare parameters this way, one hopes that NSAS will be recovered in all other 
correlation functions, thus reproducing for example the results of current algebra [Q. 

This situation is unsatisfactory for several reasons. On the theoretical level, one 
would like to have a true non-perturbative proof of the restoration of NSAS. More- 
over, the construction of a lattice model of QCD where the existence of the chiral 
limit can (a) be proved and (b) does not require any fine tuning, is important for 
practical reasons. To elucidate the importance of such framework, we can mention 
for example the problems involved in measuring weak matrix elements on the lat- 
tice with Wilson fermions p, 0. The fine tuning problem is not over when the bare 
masses have been fixed. Because of the hard breaking of the axial symmetries, the 
definition of renormalized four fermion operators which are necessary for the compu- 
tation of weak decays, involves additional fine-tuning. Moreover, some of the relevant 
operator mixings receive genuinely non-perturbative contributions, which cannot be 
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determined even in principle by short distance expansions such as the OPE. 

An alternative lattice formulation which does enjoy a certain degree of axial 
symmetry is the staggered fermions formulation. The staggered fermions action has 
a Uv{l) X f/A(l) symmetry, and the exact Ua{1) can be used for a better determination 
of meson decay constants . However, this formulation also has several drawbacks. 
The number of flavours must be equal to four, and disentangling spacetime and flavour 
symmetries is a non-trivial issue. Moreover, in the continuum limit one expect to 
recover the full SU{4) axial flavour symmetry, but on the lattice only one of the 
corresponding currents is conserved. The other fourteen currents presumably suffer 
from the same problems as in the case of Wilson fermions. 

In this paper we present a new lattice formulation of QCD with a very mild break- 
ing of all non-singlet axial symmetries. The formulation is based on the introduction 
of many (in the chiral limit infinitely many) heavy "regulator" fields [|lT], 0, |l^ . More 



specifically, we use a variant |T^ of Kaplan's proposal |T^ to realize light ordinary 
fermions as zero modes bound to four dimensional defects in a theory of massive five 
dimensional Dirac fermions [12-21]. The right-handed (RH) and left-handed (LH) 
components of the physical quark arise as surface states on opposite boundaries of a 
five dimensional slab with free boundary conditions in the fifth direction. One five 
dimensional fermion field is needed for every physical quark. The surface fermions 
scheme has also been discussed recently in ref. 

While Kaplan intended to propose a solution to the long-standing problem of 
defining chiral gauge theories on the lattice, the feasibility of reaching this goal is still 
unclear [16-19,22]. But the advantages of using chiral defect fermions for lattice QCD 
are obvious. In particular, one can easily show that, in the limit where the width of 
the five dimensional slab tends to infinity, quark masses undergo only multiplicative 
renormalization to all orders in perturbation theory pi|. 

The reason for the absence of additive 0{l/a) corrections to quark masses, is the 
vanishing of the overlap between the RH and LH components of the quark's wave 
function. At tree level, the tail of the RH wave function goes like 

{1-May. (1.1) 

Here s is the fifth coordinate which takes the values s = 1, . . . , 2A^. The Dirac mass 
M that appears in the five dimensional fermion action obeys < Ma < 1. A similar 



expression applies to the LH component, but with s in eq. ( |1.1| ) replaced by 2A^ — s. 
Thus, the perturbative overlap of the LH and RH components vanishes exponentially 
with increasing A^. 

One can also introduce direct couplings between the LH and RH components 
of every quark, by adding links that couple the layers s = 1 and s = 2N. These 
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couplings are controlled by dimensionful parameters m^, where z is a flavour index. 
As shown in ref. , the rrii play the role of multiplicatively renormalized quark 
masses. 

Taking advantage of the special properties of the model, one can define axial 
currents whose divergences, for rrii = 0, are completely localized on the two middle 
layers s = N and s = N + 1. As a result, the anomalous term in the Ward identities 
of NSAS is governed by the small tail of the quark's wave function at the center 
of the five dimensional slab. On the other hand, the divergence of the singlet axial 
current can couple to two gluons, giving rise in the limit — > oo to the expected 



axial anomaly [^, [Ty, |15 



In this paper we extend the investigation of axial properties to include also non- 
perturbative effects. Our main result is that NSAS become exact in the limit A^ oo. 
To our knowledge, this is the first non-perturbative proof of the restoration of NSAS. 

The restoration of NSAS occurs before the continuum limit is taken. The limiting 
'W = oo" formalism should be regarded as non-perturbative regularization of QCD 
which is maximally symmetric under axial transformations. Since the 'W = oo" 
formulation preserves NSAS while reproducing correctly the singlet anomaly, it cannot 
be a strictly local regulator. The non-locality arises from integrating out infinitely 
many heavy four- dimensional fields, and we believe that it is mild enough not to 
jeopardize the consistency of the continuum limit. 

This paper is organized as follows. In sect. 2 we give the definition of the model. 
Apart from the obvious gauge and fermion fields, the model includes a set of massive 
five dimensional scalar fields. These fields are necessary to cancel out some lattice 



artifacts of the five dimensional fermions [0, |T^. Such scalar fields have often been 
called Pauli-Villars (PV) fields in the literature, and we will continue to use this 
terminology here. But it should be stressed that the scalar action is non-negative, 
and so the partition function of the PV fields is well-defined. The peculiar property 
of the PV fields is that, apart from a different choice of boundary conditions, their 
action is the square of the fermionic action. 

In sect. 3 we develop a transfer matrix formalism to represent the fermionic 
partition function as well as correlation functions. The transfer matrix technique was 
introduced in the context of chiral defect fermions by Narayanan and Neuberger |jl6| 



and it relies on the work of Liischer It proves particularly convenient for the 

investigation of the model. Formulae are given both for finite A^ and for the limiting 
case N ^ oo. 

Sect. 4 contains a discussion of the effective action S^{U) obtained by integrating 
out the fermion and PV fields and taking the limit N ^ oo. An interesting result 
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is that, while being always real, exp{—S^{U)} is not necessarily positive. This 
behaviour can be explained on the basis of familiar instanton results. 

The results of sect. 3 and 4 are valid for a fixed background gauge field, provided 
the hamiltonian Ti = — log T where T is the transfer matrix has a unique ground 
state. In sect. 5 we derive an analytical criterion eq. (|5.3|) which selects those gauge 
field configurations where Ti has an exact zero mode and, hence, the ground state 
is degenerate. Such backgrounds allow for unsuppressed fermionic (and PV) propa- 
gation across in the fifth dimension, and they contribute to the anomalous term in 
non-singlet axial Ward identities. The possibility of having light states inside the 
five dimensional balk is related to the unconventional relative sign of the mass and 
Wilson terms in the fermion action. (The hopping parameter is supercritical ). A 
simple example of a configuration which supports anomalously light balk states is the 
dynamical domain wall. 

Sect. 6 contains the main results of this paper. Let Z{gQ, L, N) be the partition 
function on a finite five dimensional lattice. The N ^ oo limit of the partition 
function is 



Both the four dimensional lattice and the bare coupling qq are kept finite. We 
prove that 



Sg{U) and S^{U) are defined in eqs. ( p^) and ( p.l9| ) respectively. In other words, 
one can interchange the order of integration over the group variables and the N oo 
limit. This result follows from the s- independence of the gauge field, compactness of 
the gauge field configuration space and eq. (|5.3| ), which implies that exact zero modes 
of 7i exist only on a zero measure subspace. 

Repeating the same steps for non-singlet axial Ward identities, we prove that the 
anomalous term in every Ward identity vanishes in the limit ^ oo, at fixed values 
of L and Qq. As expected, the same analysis gives rise to a non-zero expression for 
the singlet anomaly. We stress that NSAS are recovered before the continuum limit 
(7o ^ is taken. Such a non-trivial result is possible because the size of the unphysical 
fifth dimension A^ is sent to infinity. This is reflected in the operator expression for 
conserved currents, which contains an infinite sum over the s-coordinate. 

A difficulty with the transfer matrix formalism, it that its efficient implementation 
in numerical simulations would require the development of new techniques. Instead, 
one may choose to put the fermions on a finite five dimensional lattice. One should 
then choose an optimal value for A^ subject to the constraints dictated by computer 



Z^{go,L)= hm Z{go,L,N). 





(1.3) 
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performance. To this end, it is important to have a reahstic estimate of the magnitude 
of anomalous effects on a finite five dimensional lattice. 

We believe that the bounds used in proving the existence of the chiral limit, 
highly overestimate the true magnitude of anomalous terms. A detailed study of the 
issue is beyond the scope of this paper. However, we have decided to give in Sect. 7 
a short heuristic discussion, which is meant to give the reader some feeling about 
the plausible magnitude of the anomalous term. As we explain, we believe that the 
magnitude of anomalous effects may turn out to be numerically very small already 
for currently accessible five dimensional lattices. Finally, some technical details are 
relegated to two appendices. 



2. Definition of the model 

In this section we give the definition of lattice QCD with the surface fermions 



variant of chiral defect fermions . We also define physical quark operators as well as 
axial and vector currents |18[ appropriate for the model. Most of the ingredients have 
been introduced previously, and we give them here to make the present exposition 
self-contained. 

For definiteness, we take the physical case of four dimensions. This means that 
the fermion and Pauli-Villars (PV) fields live on five dimensional lattices, whereas the 
gluon fields are four dimensional. The ordinary four coordinates, labeled x^, range 
from 1 to L, whereas the extra coordinate takes the values s = 1, . . . ,2N for the 
fermionic lattice. The PV lattice is only half as big, with s ranging from 1 to A^. The 
preferred boundary conditions in the four ordinary dimensions are periodic or anti- 
periodic. Free boundary condition in these directions would result in extra unwanted 
light states which can propagate along the spatial boundaries. 

The above scheme is realized by requiring that the link variables in the fermion 
and PV action obey Ux,s,5 = 1 and Ux,s,^i = Ux,^ independently of s. The topology 
of the fifth dimension is taken to be a circle, but the couplings which reside on the 
links connecting the layers s = 2N and s = 1 are proportional to a parameter — mj 
(i = 1, . . . , Nf is a flavour index. Also, we henceforth set the lattice spacing to a = 1). 
The case rrii = 1 corresponds to antiperiodic boundary conditions, where the model 
supports no light fermionic state. The case = corresponds to open boundaries, 
and it should give rise to the physics of QCD with massless quarks by taking first the 
limit N ^ oo and then the continuum limit. 
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The partition function is 



Z{go,L,N, rrii) 

2N „ N 



= n n / ^^^-,^11 / di'.^sd^.^s n / (2-i) 

X \ ^l 3=1'' s'=i J 

The action is given by: 

5 = Sg{U) + Sf^, tfj, U) + Spvi'pK 0, U) . (2.2) 

Colour, flavour and Dirac indices will be suppressed throughout this paper unless 
explicitly stated otherwise. We remind the reader that the PV fields carry the same 
set of indices as the fermion fields. Also, we will usually write T>U = JJ dU etc. as a 
shorthand for the corresponding measure. 

SciU) is the pure gauge part of the action. The results of this paper generalize 
to any compact Lie group with a non-negative lattice action, which reduces to the 
appropriate gauge field action in the classical continuum limit. To avoid irrelevant 
notational complications, we will assume that the gauge group is G = SU{Nc). In 
eq. (^?T|) , dU is the normahzed, invariant group measure, and Sg{U) is the the usual 
sum over plaquettes 

Sg = \Y1 E Retr(/-f/,,^,). (2.3) 

Here Ux,^u is the plaquette variable. 

The fermion and PV actions contain a sum over all flavours. The only difference 
between various flavours can be in the mass parameter rrii. We give below the one 
flavour action. The fermionic part has the following form 

Sj.{^, ^,U) = - i^xADF)x,s;y,s' i^y,s' , (2.4) 

x,y,s,s' 

where the fermionic matrix is deflned by 

{DF)x,s;y,s' = ^s,s'DI + Sx^yDj-,,, , (2.5) 



^iy = ^E((l+7/.)t^x,A+/i,y + (l-7M)f^yVx-A,yy 

+ (M-4)5,,,, (2.6) 



Pr 52,s' - mPi 62N,s' - Sl^s' , 5 = 1, 

Dts' = { Pr 5s+i,s' + Pl 6s-i,s' - 5s,s' , K s < 2iV , (2.7) 

-ttiPr 6i^s' + Pl S2n~i,s' - S2n,s' , s = 2N . 
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and 

Pr,l = 1{1±1,). (2.8) 

Notice that D;^^, is independent of the gauge field. Also, apart from the uncon- 
ventional sign of the mass term, DI y is the usual four dimensional gauge covariant 
Dirac operator for massive Wilson fermions. 

When rrii = the spectrum of surface states contains one RH Weyl fermion near 
the boundary s = 1 and one LH Weyl fermion near the other boundary for every five 
dimensional fermion field. These Weyl fermions have the same coupling to the gauge 
field, and so they in fact describe Nf quarks. If we ignore the exponentially small 
overlap between the tails of the LH and RH surface states, then these states describe 
massless quarks. Switching mj on, we now mix the RH and LH components of each 
quark and provide it with a Dirac mass which is proportional to mj. At tree level 
one has 

mi = M(2 - M)mi . (2.9) 



The PV fields |Tl], [T^ are needed to cancel the contribution of heavy fermion 
modes to the effective action Scs{U). This contribution, while formally being local 
in the continuum limit, is proportional to A^. (Every five dimensional fermion field 
describes one light quark field and 2A^ — 1 four dimensional fields whose mass is of the 
order of the cutoff). If one does not subtract the contribution of the massive fields 
by hand, that lattice artifact will dominate the effective action in the limit N oo. 

Let us denote the dependence of the Dirac operator in eq. ( ^.5] ) on rrii and on the 
number of sites in the s-direction by Dp = Dp{2N,mi). The PV fields live on a five 
dimensional lattice with A^ sites in the s-direction, and using the above notation, the 
PV action is 

Spvi(j)\ <P,U)= ^is DUN. l).,s;z,s" Dp{N, l).,s";y,s' <Py,s' ■ (2.10) 



The second order operator in eq. ( p. 101 ) is the square of the Dirac operator on a 



smaller lattice. The choice rrii = 1 for the PV fields prevents the appearance of light 
scalar modes on the layers s = 1 and s = N. As will be shown below, this choice of 
Spv ensures that exp{— Scfr} remains finite in the limit N —>■ oo. 

Let us denote by 71 the refiection relative to the hyper-plane s = A^ + 1/2 

7llpx,s = 1px,2N+l-s ■ (2.11) 

The Dirac operator Dp of eq. ( P3| ) satisfies the following identity 

757^DF757^ = I^^. (2.12) 



This identity is a generalization of a similar relation obeyed by the four dimensional 
Dirac operator for Wilson fermion, which reads 

75^1175 = (2.13) 

As in the case of Wilson fermions, eq. ( p.l2| ) plays an important role in establishing 
the positivity of pion correlators (see App. B). 

Eq. (|2.12| ) implies that the operator 'j^TZDp is hermitian. One has det {'y^TZ) = 1 
trivially, and so 

det (Dp) = det {-i^UDp) • (2.14) 

As a result, the fermionic determinant is real. However, one cannot conclude that the 
fermionic determinant is necessarily positive. The fermionic determinant can in fact 
change sign, due to level crossing of an odd number of states (see below). 

We comment in passing that one can use the hermitian operator '-y^TZDp in the 
definition of the fermionic action instead oi Dp. This is facilitated by the unitary 
change of variables (recall that if) and if) are independent variables in euclidean space) 

^ ^ ^' = ^757^. (2.15) 



A similar definition (involving 75 only) has been found useful in the construction of 
the euclidean partition function of continuum supersymmetric QCD |^^. Notice that 
the uneven treatment of ip and ip in eq. (|2.15|) gives rise to unconventional expressions 
for the axial and vector currents. Since we will not use this formulation any further 
we do not give the details here. 

In the rest of this paper we assume rrii = m, i = 1, . . . , Nf. The five dimen- 
sional fermion action is invariant under global U{Nf) symmetry. The conserved five 
dimensional current has the following components. For /x = 1, . . . , 4 

f^ix,s) = i(t^^,,(l + 7M)f^x,MA>x+/i,. 

-V^.+A,^(l - 7/.)t/.VA>.,.) , l<s<2N. (2.16) 
As for the fifth component, we define 

^ ' \ 1px,2NPR>^'''ipx,l-i^x,lPL>^'''ipx,2N , S = 2N . 

This five dimensional current satisfies the continuity equation 

-j^ix,l)-rnjl{x,2N), s = l, 
J2^,r,{x,s) = { -A,jl{x,s), l<s<2N, (2.18) 

j^{x, 2N -1) + mj^{x, 2N), s = 2N . 
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Here 

A^/(x, s) = fix, s) - fix - /i, s) , (2.19) 

A,fix,s) = fix,s)-fix,s-l). (2.20) 

A" is a flavour symmetry generator. Notice the special form of the boundary terms 
in the continuity equation. 

We now give the definitions of four dimensional vector and axial currents [p!8| . 
There is a unique set of conserved vector currents, given by 

2N 

V,%x)=Y.j;ix,s). (2.21) 



Conservation of the vector current follows from eq. (|2.18| ). 

On the other hand, there is a lot of arbitrariness in defining axial transformations 
in the model. Any transformation which assigns opposite charges to the LH and RH 
chiral modes will reduce to the familiar axial transformation in the continuum limit. 

Here we take advantage of the global separation of the LH and RH modes in 
the s-direction. We define our axial transformation to act vectorially on a given four 
dimensional layer, but we assign opposite charges to fermions in the two half-spaces 

= +^g(s)A>,,,, (2.22) 
Sl'4^.,s = -iqis)^^X, (2.23) 



where 



qis) = { ' - - ' 2.24 
^ ' ^ -I, N < s<2N . ^ ' 



The corresponding axial currents are 

^» = - E sign(iV - s + -) j;(a:, s) . (2.25) 

For m = 0, the non-invariance of the action under the transformation ( p.22|) 
resides entirely in the coupling between the layers s = N and s = + 1. For m ^ 0, 
there is an additional contribution coming from the direct coupling between the layers 
s = 1 and s = 2N. As a result, the axial currents satisfy the following divergence 
equations 

A,A^^ix) = 2mJ^ix) + 2J^^ix), (2.26) 

where 

J^ix) = f,ix,2N) , (2.27) 
J^^ix)=r,ix,N). (2.28) 
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In order to understand the physical content of eq. ( p.26 ) let us define quark 
operators as follows 



Qx = PriPx,i + PliPx,2N 

Qx = ^x,2NPR + i^x,lPL. (2.29) 

These operators have a finite overlap with the surface states in the chiral limit —>■ 
oo, and so they can play the role of bare quark fields. There are of course many other 
operators which are localized near the boundaries and, hence, can interpolate quark 
states. Our choice eq. ( |2.29|) is the simplest possible one, and it leads to considerable 



simplification of the expressions for correlation functions. 
In terms of the quark fields, J5 takes the familiar form 

J5"(x) =g,75A'^g,. (2.30) 

Thus, up to the finite normalization factor in eq. ( p.9| ), the J5 term on the r.h.s. 



of eq. (|2.26|) is the expected contribution from a classical mass term. J^^ represents 
an additional, quantum breaking term. Below we will be interested in axial Ward 
identities of the general form 



+2{j^x)0{y,,y2,..)) 

0(1/1, 2/2,...)), (2.31) 

where is a non-singlet axial current. Our aim is to show that, if 0{yi,y2, ■ ■ ■) 
involves only the quark operators of eq. (|2.29|) then the anomalous term 



J^^ix)0iy,,y2,...)) (2.32) 

vanishes in the limit N ^ 00. An especially important case is when the operator O 
is itself a pseudo-scalar density O = J^iy). The corresponding Ward identity 

/\,{AI{x) Jl{y)) = 2m{.r,{x) Jl{y)) 

+2{r,^{x) Jl{y) 



governs the pion mass. 



6,^y{qy{\\\'}qy), (2.33) 
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3. Transfer matrix formalism 



It is convenient to discuss non-perturbative effects using the transfer matrix for- 
malism. This technique was adapted to domain wall fermions in ref. [jl6|. Here we 



in fact have a constant five dimensional mass M, and the only deviation compared 
to ref. ||23| is in the mj-dependent value of the couplings on the links connecting the 
layers s = 1 and s = 2N. In this section and the next one we will assume that the 
background gauge field is fixed. 



A simple generalization of the result of ref. gives rise to the following second 



quantized expression for the Grassmann path integral 

det DF{2N,m) = J V'^Vifj e-^^^'^'^'^^ 

= (detS)2^trT2^0(m). (3.1) 

The second quantized transfer matrix T acts in a Fock space spanned by the 
action of fermionic creation operators al. on a vacuum state |0) annihilated by cix- 
The operators and aj. satisfy canonical anticommutation relations. They live on 
the sites of a four dimensional lattice of size L^, and they also carry Dirac, colour 
and flavour indices which we have suppressed. The decomposition of into RH and 
LH components is 

c 



The transfer matrix is defined by 



where 



T = e-^^^^ , (3.3) 



(3.2) 



^ ~ \ C+B-^ C+B-^C + B ) ^ ' ^ 

Bx,y = (5 - M)6x,y -^Yl [^X+fl,yUx,fM + Sx-f,,yUl^ , (3.5) 
^x,y = 9 X! [^x+il,yUx,,_i - Sx-fi.,yU^^^ CTfj, , (3.6) 

and = {i, a). An important identity is 

e"^ = KK^ , (3.7) 

where 

( B-^l^ , 
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The last equation implies that the matrix operator H is well-defined. Notice also that 
Z^ll can be expressed in terms of B and C as follows 




D'=\ . . • (3.9) 



The operator 0{m) contains all the m-dependence, and it is given by 

0{m) = Yiicncl + mcicn){dndi + mdldn) . (3.10) 

n 

In this equation, n is a generic name for all indices. The special cases m = and m = 1 
deserve special attention. For m = 1, 0{m) becomes the identity operator, whereas 
for m = it is a projection operator on a different ground state |0') annihilated by 
all the c-s and d-s. The relation between |0) and |0') is 

o>=n4|o'). (3.11) 

n 

Both |0) and |0') are "kinematical" ground states which can serve as convenient 
reference states in the construction of the Fock space. But none of them play a special 
role dynamically. Since we use creation and annihilation operators in the coordinate 
basis, both |0) and |0') are very different from the filled Dirac sea even in the case of 
free fermions. 

Like D", the hermitian operator H is a Nt x Nt matrix, where Nt = ANcNfL^. 
Let R be the unitary matrix which diagonalizes H 

HmlRln = EnRmn ■ (3-12) 



Corresponding to the splitting into two-by-two matrices in Dirac space in eq. 
we write R as two ^Nt x A^^ matrices P and Q. We will assume that the columns 
of R are ordered according to their eigenvalues, with the negative eigenvalues on the 
left. Accordingly, we further decomposed P and Q into submatrices which contain 
the positive and negative eigenvectors 



R=\ ^ \ ■ (3.13) 




P'F and Q"^ are ^Nt x matrices, where + A^" = A^^. 

The ground state of of the second quantized operator Ti. = a)Hd is obtained by 
filling all negative energy states 

N- 

1=1 
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In the limit — > oo, T is proportional to a projector on the ground state of Ti. 



T 



2N 



N ^oo 



(3.15) 



where 



N- 



Amax = exp | - ^ ^ . (3.16) 
We now turn to the scalar PV fields. The action is bilinear in these fields, and so 



(deti?)-2^ (trT^) 



-2 



(3.17) 



In going from the first to the second row we have used eq. ( p.l|) and substituted 
m = 1. 

The effective action Scs{U) is defined by integrating out both fermion and PV 
fields. Using eqs. (|0| ) and (|3.17|) we find 



In the limit N ^ oo the effective action becomes 



(3.18) 



exp{-5~} 



lim expl-^cflf} 

Af— »oo 



H 



0{m] 



H 



(3.19) 



Eq. ( p. 191 ) is valid for a fixed background field provided the ground state of H 
is non-degenerate. The limiting expression eq. ( |3.19| ) is completely well defined here, 
and it is free of any subtleties of the kind encountered in trying to define chiral gauge 
theories on the lattice using chiral defect fermions ]TB|. 



4. The effective action 

As a first step, we want to study the behaviour of {Oh \ 0{m) \0h) under various 
conditions. The physically interesting case is m <^ 1. In the case m = 1 one has 
{Oh\ \0h) = 1- The reason for this trivial result is the subtraction of the balk 
effect through the PV fields. By contrast, in the opposite limit m = one has 



exp{-S^} 



H 



0' 



m = . 



(4.1) 



Thus, in a model of massless chiral defect fermions, the physical information is not 
in the trace of the transfer matrix but, rather, in the overlap of its ground state with 
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some other state. This was first found by Narayanan and Neuberger for domain wall 
fermions, where the overlap formula reads WE 



exp{ (domain wall)} = (Oh+ 0//_) . (4.2) 

Here H± are the hamiltonians corresponding to the two sides of the domain wall. 
Hence, one has to compare two dynamical ground states. But all the non-trivial 
dynamics of the model is contained in the if+ hamiltonian. In the continuum limit, 
eqs. (|4.1| ) and ( [4.2|) should describe the same physics. It is therefore advantageous to 
work with the surface fermion scheme, where one has to calculate the overlap of \0h) 
with a fixed reference state |0'). 

An explicit expression for the m = overlap eq. ([4. 1|) can be easily written 
down |T6[. We first comment that for free fermions, as well as for perturbative 



gauge field configurations, the numbers of positive and negative eigenvalues of H 
are equal A^^ = jNf. In this case P"^ and are square matrices. Using eq. ( p.ll|) 
it follows that only the cl-dependent terms in eq. ( ^.14[ ) contribute to the overlap 



eq. ([4.1| ). Taking into account the anticommuting character of the fermionic operators 
one arrives at 



H 



0') = det Q- . (4.3) 



The phase ambiguity in defining the columns of is irrelevant, because only the 
modulus of detQ" enters eq. (|4.1|). 

For non-perturbative configurations there may be level crossing, resulting in 
A^^ 7^ |Af. In this case the m = overlap vanishes identically. As discussed in 
ref. [|l^ this is a welcomed phenomenon, which signals that the chiral defect fermion 
model can reproduce instanton effects. 

We now want to generalize the explicit expression for {Oh \ 0{m) \0h) to m 7^ 0. 
We first notice that 0{m) can be expanded as 



0{m) = Y.m' E IT-A ■ ■ ■ 44. • • • ■ ■ ■ • • • 4 • (4-4) 

A:=0 l+n=k 

In eq. ([4.4|) , multiplies a sum over orthogonal projection operators whose number 
grows like N^/k\. Recall that Nt is proportional to the four- volume V = L^. In 
calculating the m-expansion of the effective action, at the k-th order we will therefore 
encounter 0{V^/k\) terms, where the magnitude of each term is bounded by m''. 
This is in agreement with the anticipated behaviour of a system which undergoes 
spontaneous symmetry breaking. The finite volume partition function is analytic in 
m, but in the infinite volume limit singularities may appear because the product mV 
diverges. 
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Returning to the calculation of (Oi^| 0{m) \0h), let us denote 

iNt - N- 



(4.5) 



The first term in the expansion eq. ( |4.4| ) which contributes to (0^1 0{m) \0h) is the 
term with k = A. This immediately implies that 



H 



0{m] 



H 



(4.6) 



Moreover, for sufficiently small m, the sign of {Oh\ 0{m) \0h) is determined by 
the sign of m^. Hence, (0^1 0{m) \0h) will be negative for m < and odd A. For 
example, in the one fiavour case {Oh \ 0{m) \0h) is negative for m < if an odd 
number of level crossing took place. 

This behaviour is not unexpected and, in fact, it is in agreement with the familiar 
instanton result p5[. In the one fiavour case, the fermionic determinant in an instan- 
ton background is proportional to m, and so it changes sign if m does. The gauge 
field's effective measure exp{— 5^ — 5*^} is therefore real, but not always positive. 
We still expect the partition function to be strictly positive when one approaches the 
continuum limit, because configurations with a non-zero topological charge are rare. 
But we are unable to prove the positivity of the partition function in a completely 
general way. Special cases where exp{— S*^} is strictly positive include m > 0, or 
m 7^ and even Nf. 

An explicit expression for {Oh \ 0(171) \0h) is more easily obtained using the defi- 
nition eq. (13.101 ). Straightforward application of the canonical anticommutation rela- 
tions gives rise to 



N- 



i=l 



(4.7) 



Introducing the Nt x N matrix 



R-(m) 




we find 



IH 



0(m) 



IH 



det R^-(1)K 



In the case m > one can use the relation 0{m) = 0{m2)0{m2) to write 



IH 



0{m) 



det R^-(m2) R- 



(4i 



(4.9) 



(4.10) 



which is manifestly positive. 
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The transfer matrix formalism can also be used to write down expressions for 
correlation functions. The correlation functions for the quark operators ( [^.29 ) take 
a particularly simple form. The rules for their construction, as well as some explicit 
examples, are given in Appendix A. The transfer matrix formula for axial Ward 
identities is discussed also in Sect. 6. 



5. The dynamical domain wall 

In the previous section we discussed some properties of the effective action under 
the assumption that the ground state of 7i is non-degenerate. The ground state will 
be degenerate if Ti has an exact zero mode or, equivalently, if the matrix has a 
unit eigenvalue 

e-^V'o = ^o. (5.1) 

If eq. ( |5.1| ) holds for a particular background field we expect strong correlations 
accross the five dimensional slab. In particular, the anomalous term in NSAS Ward 
identities receives its dominant contribution from such configurations. It is therefore 
important to identify these configurations and to understand their properties. 

The matrix operator is non-local, and to understand better the physical 



content of eq. (|0| ) we want to find a simpler equation satisfied by the zero mode ifjQ. 



Using eqs. (fj.7D, ( |3.8D and ( ^.91 ) it is straightforward to show that 

= {K^-K-^)%l)o 

= fi-i/S^"^o. (5.2) 

Since 5 is a positive definite operator, we conclude that ipo is a zero mode of the 
hermitian operator 75-D". Therefore, a necessary and sufficient condition for the 
existence of a zero mode is 

det(75L'll) = 0. (5.3) 
Notice that the l.h.s. of eq. ( ^.3|) is a polynomial in the group variables Ux,^- 



Since is a massive lattice Dirac operator, one may ask whether eq. ( |5.1| ) has 
any non-trivial solutions at all. An indirect way to argue that such solutions should 
exist, is to observe that this is a necessary condition for (Oj^l 0{m) \Qh) to reproduce 
the instanton results. The vanishing of this expectation value for m = requires that 
level crossing should occur in the spectrum of H . At the crossing point one has a 



solution of eq. ( ^.1|) . This observation was made by Narayanan and Neuberger p6 
who also found numerically solutions of eq. ( |5.1| ). 
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The reason why solutions to eq. (|5.1|) exist in spite of the mass term present 
in D'I = DII(M), is the unconventional sign of that mass term. For comparison, the 
conventional Dirac operator for a massive Wilson fermion is D" (— M) in our notation. 
The "wrong" sign of the mass term relative to the Wilson term, implies that the 
sum of the two terms is not a positive definite operator. As a result, Vafa-Witten 
bounds p6[ cannot be established here. On the other hand, if one were to choose 



the conventional relative sign, then a Vafa-Witten bound could be establish for the 
propagation of fermions in all directions. This, in turn, would imply the absence of 
any light states in the model. Indeed, one can easily check that the massless surface 
modes disappear for M < 0. We comment that in the domain wall case, too, all 
the non-trivial dynamics occurs on that side of the wall where one has the "wrong" 
relative sign. 

In the absence of gauge fields, as well as in perturbation theory, the only zero 
modes of the five dimensional Dirac operator Dp are the ones discussed in refs. |]12 



21] , |TJ]. But with dynamical gauge fields, other zero modes can exist. Since the gauge 
field is s-independent, any zero mode of the four dimensional operator Z^H will also be 
a zero mode of the five dimensional operator Dp. This is true up to boundary effects, 
which can be ignored here because the new type of zero modes in not localized in the 
s-direction. 

An interesting example is provided by the dynamical domain wall. For simplicity 
we consider here a U{1) gauge group. (The same configuration exists also for SU{2). 
For general SU{N) one can simply embed the below configuration in some SU(2) 
subgroup). We consider the following configuration of link variables. We let Ux,fj. = 1 
for fi = 1,3, 4. For Ux,2 we take 

U^,2={ I' ^^^^1' (5.4) 

-1 , Xi > . 



Notice that this is essentially a two dimensional configuration. Eq. ( p.4|) describes a 
wall of magnetic flux, with one unit of flux going through each plaquette to the left 
of the hyperplane xi = x?. 

We now consider the ansatz 

^o = ^(l-7i)*o(a;i). (5.5) 

One can easily check that the following is a zero mode of Z?" 

,T, r 1 / (l-A!f)'°-", x,<xl 
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This is the simplest case of a dynamically generated zero mode. Other gauge 
field configurations with a topological character may also support zero modes of D". 
These include for example lattice analogs of the static string-like singularity discussed 



Finally, we recall the close relationship between the fermionic and PV action. 
Whenever the hermitian operator 75-D" has light states, the same will be true for 
the PV fields. Therefore, anomalously strong correlations in the s-direction arise 
simultaneously for the fermions and the PV fields. However, the two contributions 
to the anomalous term do not cancel each other in general. 

In the Ward identity eq. (|2.33|) which governs the pion mass, the anomalous term 
is the correlator of two pseudoscalar densities. As a generalization of similar results 
for Wilson fermions, one can prove in the present model the positivity of a whole 
family of correlators of pseudoscalar densities. The details can be found in App. B. 
As a result, the anomalous term in eq. ( |2.33| ) is strictly positive for finite values of 
N and qq. The vanishing of the anomalous term in the limit N 00 can only arise 
from the vanishing of its absolute value on almost the entire gauge field configuration 
space. 

6. The chiral limit 

The characterization (|5.3| ) of gauge field configurations which support exact fer- 
mionic zero modes, allows us to prove the following important result (see eq. (p..2[)) 



In technical terms, eq. (|6.1|) means that we can interchange the order of the N 00 
limit and the integration over the group variables. As we discuss below, the physical 
content of eq. ( |6.1| ) is that correlation functions obey clustering in the s-direction. 

Let us first introduce some terminology. The gauge field configuration space is 
g = {SU{Nc))^^ . An element of G, i.e. a particular configuration of link variables 
will be denoted hj U = {f/x,/i}- We let Go <Z Q he the subspace of all configurations 
which satisfy condition (|5.3|) . 

In this section we make the technical assumption < m < 1. (The reader should 
not confuse this condition with the similarly looking one < M < 1, which we 
enforce for entirely different reasons). The condition m > ensures the positivity 
of {Oh\ 0{m) \0h)- In the continuum limit we expect to recover the symmetry under 
m — m, and so choosing a non- negative m should not lead to any restrictions on 
the physical content of the model. The condition |m| < 1 implies || 0{m) || < 1. This 



in ref. ||2^. The zero modes observed in ref. are actually of that type. 




(6.1) 
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gives rise to less cumbersome expressions for some of the bounds below. As we have 
already explained, the physically interesting case is |m| <^ 1. 

The proof of eq. ( |6.1|) is simple, and it relies on the following ingredients, (a) The 
gauge field is s-independent, and the configuration space Q is compact, (b) According 
to standard results in calculus, the subspace Qq has a zero measure. 

Consider an element U G Q — Go, and let En be the eigenvalues of H{U). We 
define 

Eo(W) =min{|E„|}. (6.2) 

Since H{U) is a finite dimensional matrix, Eq{U) is well-defined, and since U G Q — Qo, 
one has Eq{U) > 0. 

We use this information to put a bound on the difference between the r.h.s. of 
eq. (|3.18| ) and the r.h.s. of eq. ( |3.19|) . Separating the ground state contribution from 
the rest we find 



tr T2^C(m)/ (trT^) - {0h\O{ 
A^t(trT^)-^-l 



H 



< 



2N 



+ 



< 2 



Nt 



(6.3a) 

(6.3b) 
(6.3c) 



The notation tr ' means that the ground state contribution is excluded. Notice that 
2^* is the dimensionality of the fermionic Fock space. The last row of the above 
inequality highly overestimates its first row. Nevertheless, it will be sufficient for our 
purpose, because it implies that expression (|6.3a| ) vanishes in the limit iV — ^ cxd for 
Eo{U) > 0. 

We now have to show that for arbitrary e > 0, there exist such that 



Z{go,L,N,m) - J VU e'^^'-^^ (0i^|O(m) 



H 



< e 



for every N > N^. To this end, we divide the group integration into two parts 



VU 



VU+ / VU 

g, Jg-g, 



(6.4) 



(6.5) 



Qe is an open covering of Qq whose volume is less than e/4. Such an open covering 
can always be found. Furthermore, the r.h.s. of both eqs. (|3.18 ) and ( 3. IS ) is always 
bounded by one. Hence, the contribution to the l.h.s. of inequality ( |6.4| ) coming from 
the integration over is bounded by e/2. 
Now let us denote 



Eo = mm{Eo{U) \ U eG-G,} 



(6.6) 
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Since ^ — is compact, the minimum exists and satisfies > 0. It now follows 
that the contribution to the l.h.s. of inequality ( p.4| ) coming from the integration over 
Q — Qe bounded by expression ( |6.3c|) with Eq{U) replaced by Eq. The existence 
of such that inequality ( |6.4|) holds now follows from that fact that Eq > Q. This 
completes the proof of eq. (|6.1| ). 

Similar results can be established for correlation functions. In order to prove 
the restoration of NSAS in the limit — ^> cxd, we have to show that the anomalous 
term ( ^.321 ) in the Ward identity eq. ( |2.31D vanishes for any operator 0{yi,y2, ■ ■ .) 
which is constructed solely out of the quark operators in eq. ( p.29[ ). 

Consider first a fixed background field U E Q — Qq. Suppressing the coordinates 
yi, 1/2, ■■ ■ we find in the limit iV — > oo 



N^oo 



lim Jt(x)0 



5q\ 



U,N 



H 



OLO{m)OR 



H 



(6.7) 



The subscript U indicates the we have given the expression for the correlator in a 
fixed background. Here 

jl{x)=clX'c,-dl\-d^. (6.8) 
The operator expressions for Ol and Or can be found using the rules of App. A. For 



finite A^, the difference between the expressions on the l.h.s. and the r.h.s. of eq. ( |6.7] ) 
obeys a bound analogous to inequality (6.3). The crucial observation is that the 
matrix element (0/f| j5(x) \Qh) is proportional to tr A'^, and so it vanishes identically 
for a non-singlet pseudo-scalar density. Hence 



' /U,N 



(6.9) 



The constant c is the product of the operator norms ||OlOj?|| and Integrating 
over the group variables and following the same reasoning as before we conclude that 



lim (J^Jx)O 



N 







(6.10) 



This proves the restoration of NSAS in the chiral limit N ^ oo. As an example, 
the explicit expression for the N ^ oo limit of the ward identity ( p. 331 ) is given in 
App. A. 

We conclude this section with two comments. For the singlet current, trJ 7^ 
0, and the matrix element {Qh\35{x) \Qh) is in general non-zero. As was shown in 



refs. |T8|, [T^, |T5| the axial anomaly is reproduced correctly in the continuum limit. 

Our second comment concerns the factorized form of the r.h.s. of eq. (|6.7|) . In 
the limit A^ — > 00, factorization of the expectation value of an operator product 
occurs whenever the limiting s-separation between two factors in the product tends 
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to infinity. This means that correlation functions obey clustering in the s-direction. 
Now, clustering is by itself a rather weak condition, which we normally expect to 
hold even if massless particles can propagate between two points. As is often the 
case with rigorous bounds, we believe that the actual damping of correlations in the 
s-direction is much stronger, and that the correlation length in the s-direction is finite 
in physical units, if not in lattice units. A heuristic discussion of the actual magnitude 
of anomalous effects for finite is given in the concluding section. 



7. Conclusions 

In this paper we discussed in detail a new formulation of lattice QCD which 
is based on the surface fermions scheme. Our main result is that non-singlet axial 
symmetries become exact in the chiral limit. The chiral limit is defined to be the 
limit of an infinite fifth direction, at fixed finite values of the bare coupling and the 
size of the four dimensional lattice. The vanishing of the anomalous term in all NSAS 
Ward identities implies in particular that non-singlet axial currents do not undergo 
any renormalization in the chiral limit, as should be the case for truely conserved 
currents. 



As was shown in previous works |I6|, |I5|, the singlet axial anomaly is 

correctly reproduced if one takes first the limit N ^ oo and then the continuum 
limit go ^ 0. For finite Qq, the limiting 'W = oo" formulation is therefore a new 
non-perturbative regulator of QCD which is maximally symmetric under axial trans- 
formations. This new regularization scheme is mildly non-local because we have 
integrated out an infinite number of heavy four-dimensional fields. We believe that 
the mild non-locality does not jeopardize the consistency of the continuum limit, but 
we have not investigated every possible aspect of this issue. 

The non-singlet currents defined in eq. (|2.25| ) are exactly conserved in the limit 



N —>■ (yo regardless of the value of qq. Thus, they will be conserved not only in the 
continuum limit, but also in the strong coupling limit. However, it is not necessarily 
true that in the strong coupling limit, these currents retain the physical significance 
of axial current. We do not rule out that for sufficiently strong coupling, some or all 
of the doubler modes reappear. The likely consequence would be that the currents of 
eq. (|2.25| ) become vectorial with respect to the new massless spectrum. In this case, 
the singlet current would become vectorial too, and it will be conserved in the strong 
coupling phase. 

Returning to the physically interesting limit (70 -C 1, the biggest advantage of the 
new formulation is that fine tuning is no longer needed in the fermion sector. For 



22 



example, the theoretical values of current masses are determined using weak coupling 
perturbation theory, and they involve only multiplicative renormalization. Likewise, 
operator mixing are restricted by the naive transformation properties under non- 
singlet axial symmetries, and meson decay constants can be inferred directly from 
the corresponding Ward identities. 

All the above properties become exact in the limit N —>■ oo. Looking forward 
to the implementation of the surface fermions scheme in numerical simulations, it is 
important to have a realistic estimate of the magnitude of anomalous effects on a 
finite five dimensional lattice. While sufficient for proving the restoration of NSAS 
in the limit, we believe that the rigorous bounds used in Sect. 6 represent a gross 
overestimation on the actual magnitude of anomalous effects. We have decide to 
include here a short heuristic discussion of this issue, mainly because we believe that 
the true picture is much more promising. A more detailed study is relegated to a 
separate publication. 

Our central observation is the following. As we noted previously, the naive con- 
tinuum limit of the lattice operator D^^ is a massive Dirac operator. This does not 
rule out the existence of (exact or approximate) zero modes of -D" for certain gauge 
field configurations. But it does indicate that the relevant configurations cannot in 
any sense be the discretized approximation of smooth continuum gauge fields. It may 
therefore be possible to prove the following conjecture: if a given lattice gauge field 
is the discretized approximation of a smooth continuum gauge field, then there is an 
0(1) gap in the spectrum of the (hermitian) operator 75Z)". We have intentionally 
omitted here a precise definition of what we mean by a "lattice approximation of a 
smooth continuum gauge field" . Basically what we have in mind is a condition which 
states that the gauge field's action density is very small everywhere. However, there 
may be other definitions that have the same physical content, and which are more 
convenient from a mathematical point of view. 

The above conjecture receives circumstantial evidence from the two known exam- 
ples of zero modes of DK These are the singular fiuxon of ref. |TB|, and the singular 
dynamical domain wall of Sect. 5. Both are characterize by an action density which 
is 0(1) on some n-dimensional subspace, where n = 2 and n = 3 respectively for a 
fiuxon and for a domain wall. Now, if we want the configuration to support an ap- 
proximate zero mode whose energy is -Eq -C 1, then its longitudinal extension should 
be at least 0{1/Eq). This implies that the total gauge field action for such config- 
urations should be bounded from below by Sq ~ C/^QqEq). Here C is some 0(1) 
constant. 

Typically, a single extended configuration supports a single (approximate) zero 
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mode, up to symmetry factors. In other words, level crossing in the spectrum of 
Ti occurs one at a time, and all other eigenstates are separated by a finite gap in 
the vicinity of the crossing point. Assuming this to be true in general, all factors of 
2-^* (which count the total number of eigenstates) in the inequalities of Sect. 6 can 
be dropped. Putting everything together, this suggests that the magnitude of the 
anomalous term for a given background should be 



To complete the estimate of the anomalous term, it is necessary to determine what 
are the most important singular configurations. If all configurations have n > 0, the 
most important effect should be the need to increase the longitudinal extension of the 
singular configuration with decreasing Eq. The dominant configurations, as well as 
the magnitude of the resulting anomalous term, can then be determined by a saddle 
point approximation. If, on the other hand, there exist point-like four dimensional 
configurations which support (approximate) zero modes of 75-D", then a more refined 
analysis would be needed, and one would have to estimate the phase space for such 
configurations as a function of g^. We believe that, either way, the correlation length 
in the s-direction should turn out to be finite. If it is found that, moreover, the 
correlation length in the s-direction remains finite in lattice units, then we may hope 
to obtain good results already on five dimensional lattices which can be simulated 
today. 

This research was supported in part by the Basic Research Foundation admin- 
istered by the Israel Academy of Sciences and Humanities, and by a grant from the 
United States - Israel Binational Science Foundation. 

A. Correlation functions in the N = 00 limit 

In this appendix we give the transfer matrix formulae for various correlation 
functions. The correlation functions for the quark operators ( |2.29| ) take a particularly 
simple form in the limit N = 00. For example, the expression for the quark condensate 
is 



The subscript U indicates that we give the expression for the correlator in a fixed 
background field. Another example is the two pion correlator 
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Here 



+ {0h 
'Oh 



Oh 
Oh 
Oh 



(A.2) 



TT 



TT 



(A.3) 



The arrows denote isospin. Notice that the pion operators are special cases of the 



pseudoscalar densities eq. (p.27|) . 

The general prescription for correlation functions of the quark operators eq. ( p.29| ) 
is the following. Considering qji i and qFt,L as Grassmann variables, one first reorder 
each product of quark operators such that qi and qi occur to the left of all q^ and 
(Jr. This step may result in a minus sign. One then translates the result into a matrix 
element of the form 



H 



■0{m) 



H 



(A.4) 



The operators to the left of 0{m) are obtained from the ordered product of g^-s and 
qL-s by the substitution 



Ql 

(IL 



(A.5) 



Similarly, on the right of 0{m) one makes the substitution 



d, 
— c . 



(A.6) 



All indices are left unchanged in this substitution. (The transition from the Grass- 
mann path integral to operator language involves a non-local transformation at an 



intermediate step [|23|, |T6|. But this non- locality cancels out in the final expression). 

Apart from quark operators, we may also be interested in the transfer matrix 
formulae for correlation functions that involve vector or axial currents. As an example, 
for finite N, the correlator on the l.h.s. of eq. ( p^.31 ) becomes 



2N-1 



Al{x)0{yuy2, 



u 



E sign(Ar - , - -) ^--^^ 



s=0 



where 



f(x,s) = T^fAx)T- 



(A.7) 



(A.^ 
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To obtain an explicit expression for j^{x), the different terms in eq. ( 2.1£ ) are trans- 
lated according to the following rules 



i^R,x i^L,y -dx T dl T"^ 

^L,x'^L,y -TcldlT-^ 

'^R,x^R,y -dxCy (A.9) 

The somewhat unexpected appearance of the transfer matrix T in these rules is due to 
our definition of exp{—H} eq. (|3.4|) . If instead we decided to use exp{—H'} = K'^K 



and the corresponding transfer matrix T', then there would be no factors of T' in 
resulting expression for On the other hand, the expressions for quark correlators 

would become more cumbersome. 

Convergence of the infinite sum on the l.h.s. of every Ward identity is guaranteed 
by the finiteness of the r.h.s. We comment that, if a slightly stronger form of clustering 
than the one proved in Sect. 6 holds, then 

oo 

Jim 0(yi, 2/2, • • 0)^ = E (OhI dLOim)dnj;ix, s) \0h) 

s=0 

- E {0H\j;{x,s)dLO{m)dH\0H) . (A.io) 

This equation is in particular valid if, as we have every reason to believe, the corre- 
lation length in the s-direction is finite. 



B. Inequalities 

Eq. ( p.l2| ) implies an analogous identity for fermion propagator, considered as a 
matrix. Writing the coordinates explicitly one has 

75^(0;, s; y, s'h, = G\y, 2N + 1 - s'; x,2N + 1 - s) . (B.l) 

Here the dagger refers only to the the suppressed internal indices. We will use this 
identity to prove the positivity of correlators of the following pseudoscalar densities 

K''{x,t) = + t)PR AXx, + 1 - t) 

A^ + 1 - t)PL AXx, A^ + t) . (B.2) 

Notice the special cases 

K^{x,0)=zJ^^{x), (B.3) 
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K''{x,N) = iJ^{x). (B.4) 
Using eq. (|B.1|) and inserting a minus sign for the closed fermion loop we now obtain 

(K'^ix, t)K\y, t')) = I DUe-^^^^ (det DFiU)ff I{U; x, t; y, t') , (B.5) 

where 

/ = ti [^PRG{y,N + t'-x,N + l-t)PRG\y,N + t'-x,N + l-t) 

+PRG{y,N +l-t';x,N + l-t)PLG\y,N + l-t';x,N + l-t) 

+Pl G{y, N + t';x,N + t)Pn G^{y, N + t';x,N + t) 

+ PLG{y, N + 1 - t';x, N + t)PLG\y, N + 1 - t';x, N + t)} . (B.6) 

For even Nf or for m > 0, the factor [det DfIU))^-^ is positive. We assume that 
one of these conditions is satisfied. It is now straightforward to prove the positivity 
of each term in eq. ( |B.6| ). Consider the second row as an example. Independently of 
the values of x and y, it has the generic form 

ti {PrAPlA^} = ti {p'rApIa^} 

= ti {{Pr A Pl){Pr A Pl)^} . (B.7) 

The last row is manifestly positive. 

Consider now the Ward identity ( |2.33| ) which determines the pion mass. Notice 



that 

n^{y) = K\x,N). (B.8) 
Using eq. (|B.3| ) the anomalous term in this Ward identity can be written as 

{K%x,0)K%x,N)) , (B.9) 

which is positive according to the previous discussion. The same reasoning proves 
the positivity of the two-pion correlator. 

The positivity of the two-pion correlator can also be established directly from the 
N = oo formula eq. ( |A.2|) . One should notice that the Fock space is a direct product 



of the Up and Down Fock spaces. Thanks to factorization of 0{m) = (9|(m) Oi{m), 
each matrix element in eq. ( |A.2| ) becomes the product of two complex conjugate 
matrix elements, one in the Up Fock space and one in the Down Fock space. 
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